Phase-space theory for dispersive detectors of superconducting qubits 



I. Serban, 1 - 2 E. Solano, 3 - 4 - 5 and F.K. Wilhelm 2 

1 Physics Department, Arnold Sommerfeld Center for Theoretical Physics, and Center for NanoScience, 
Ludwig-Maximilians-Universitat, Theresienstrasse 37, 80333 Munich, Germany 
S IQC and Dept. of Physics and Astronomy, University of Waterloo, 
200 University Ave W, Waterloo, ON, N2L 3G1, Canada 
3 Physics Department, Arnold Sommerfeld Center for Theoretical Physics, and Center for NanoScience, 

Ludwig-Maximilans-Universitat, Theresienstrasse 37, 80333 Munich, Germany 
* Max- Planck Institute for Quantum Optics, Hans-Kopfermann-Strasse 1, D-85748 Garching, Germany 
1 ' Seccion Fisica, Departamento de Ciencias, Pontificia Universidad Catolica del Peru, Apartado 1761, Lima, Peru 

(Dated: February 6, 2008) 

Motivated by recent experiments, we study the dynamics of a qubit quadratically coupled to 
its detector, a damped harmonic oscillator. We use a complex-environment approach, explicitly 
describing the dynamics of the qubit and the oscillator by means of their full Floquet state master 
equations in phase-space. We investigate the backaction of the environment on the measured qubit 
and explore several measurement protocols, which include a long-term full read-out cycle as well as 
schemes based on short time transfer of information between qubit and oscillator. We also show 
that the pointer becomes measurable before all information in the qubit has been lost. 
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I. INTRODUCTION 

The quantum measurement postulate is one of the 
most intriguing and historically controversial pieces of 
quantum mechanics. It usually appears as a separate 
postulate, as it introduces a non-unitary time evolution. 

On the other hand, at least in principle, qubit and de- 
tector can be described by a coupled manybody Hamil- 
tonian and thus the measurement process can be investi- 
gated using the established tools of quantum mechanics 
of open systems. Even though this does not lead to a 
solution of the fundamental measurement paradox, such 
research gives insight into the physics of quantum mea- 
surement 0, 0, i |, 0, @] . 

This basic question has also gained practical relevance 
and has become a field of experimental physics in the 
context of quantum computing. Specifically, supercon- 
ducting qubits have been proposed as building block of 
a scalable quantum computer 0, H, H, [HI • In these sys- 
tems, the detector is based on the same technology - 
small, underdamped Josephson junctions — as the de- 
vice whose state is to be detected. Thus these circuits 
are an ideal test-bed to investigate the physics of quan- 
tum measurement. Implementing a measurement which 
is fast and reliable, with a high (single-shot) resolution 
and high visibility is a topic of central importance to the 
practical implementation of these devices. 

The basic textbook version of a quantum measurement 
is based on von Neumann's postulate [ill, E3] • The state 
of the system is projected onto the eigenstatc of the ob- 
servable being measured corresponding to the eigenvalue 
being observed. This is not the only possible quantum 
measurement and has been generalized to the idea of a 
positive operator- valued measure [l3| . 

From the microscopic, Hamiltonian-based perspective, 
intensive research has been done on the measurement of 



small signals, which originated in the theory of gravita- 
tional wave detection [l4|. The main challenge has been 
to identify how signals below the limitations of the uncer- 
tainty relation of the detector can be measured - a regime 
in which the detector response is also strictly linear. This 
work has resulted in the notion of a quantum nondemo- 
lition (QND) measurement [ijj], which is the closest to a 
microscopic formulation of a von Neumann measurement. 
This result has been generalized to many other systems, 
prominently atomic physics, and also found its way to 
the superconducting qubits literature. Here, the anal- 
ogy of a tiny signal is the limit of weak coupling between 
qubit and detector. Another body of work 0, Q takes a 
more general starting point and discusses the relevance 
of pointer state and environment induced superselection. 

The measurement techniques used in superconduct- 
ing qubits are covering many of the mentioned situa- 
tions. Weak measurements can be performed using single 
electron transistors. Based on weak measurement the- 
ory, this is well understood (see Ref. 0] and the refer- 
ences therein) but only of limited use for superconducting 
qubits. There the measurements are far from projective, 
their resolution is in practice rather limited and the whole 
process is very slow. In the case of qubit, the task is not 
to amplify an arbitrarily weak signals, but to discrimi- 
nate two states in the best possible way. If the detector 
can be decoupled from the qubit when no measurement 
needs to be performed, this discrimination may involve 
strong qubit-detector coupling [la ]. 

An opposite, generic approach is to perform a switch- 
ing measurement — the detector switches out of a 
metastable state depending on the state of qubit. Switch- 
ing is a highly nonlinear phenomenon, so this type of de- 
tection is far from the weak measurement scenario. In 
most of the early generic setup, this process is a switch- 
ing of a superconducting device, e.g. a superconducting 
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quantum interference device (SQUID), from the super- 
conducting to the dissipative state [TH, E3, EH EH ■ This 
technique goes a long way, and some experiments have 
proven that the switching type of readout can achieve 
high contrast [H, [2l|. It has the drawback that it is 
not a projective readout and during the switching pro- 
cess hot quasiparticles with a long relaxation time are 
created. This limits the time between the consecutive 
measurements. Parts of this technique are well under- 
stood, such as the switching histogram (22|, the pre- 
measurement backaction (23| . and the influence of the 
shunt impedance to the SQUID [24], [2f| , but there is no 
full and single theory of this process on the same level of 
detail as the weak measurement theory. 

Recent developments of detection schemes have lead 
to vast improvements based on two innovations: instead 
of directly measuring a certain observable pertaining to 
a qubit state, one uses a pointer system, and measures 
one of its observables influenced by the state of the qubit. 
The observation is usually materialized in the frequency 
shift of an appropriate resonator, whose response to an 
external excitation links to the measurement outcome 
[IS H3, HI, HI] • These measurements offer good sensitiv- 
ity, high visibility [3(j, and fast repetition rates. They 
also allow to keep the qubit at a well-defined opera- 
tion point, although not always the optimum one. In 
many cases, the resonators in use are nonlinear — based 
on Josephson junctions. Thus, at stronger excitation, 
generic nonlinear effects can be exploited. These nonlin- 
ear effects go up to switching, which in contrast to the 
critical current switching is between two dissipationlcss 
states [3l| . Due to this performance and versatility, these 
devices also offer an ideal example for investigating the 
crossover between weak and strong measurements and 
the role of nonlinearity. 

Analyzing the properties of quantum measurement is 
an application of open quantum systems theory: the 
backaction contains a variant of projection which can be 
viewed in an ensemble as dephasing. The resolution is 
determined by the behavior of the detector under the in- 
fluence of the qubit viewed as an environment. For open 
quantum systems, a number of tools have been devel- 
oped. Most of them, prominently Born and Born-Markov 
master equations (see e.g. [32j and references therein for 
a recent review) assume weak coupling between qubit and 
environment and are hence a priori unsuitable for study- 
ing strong qubit-detector coupling. Tools for stronger 
coupling have been developed [HI, [13| but are largely 
restricted to harmonic oscillator baths and hence unsuit- 
able to treat the generally nonlinear physics of the sys- 
tems of interest. The Lindblad equation [35| is claimed to 
be valid up to strong coupling, however, due to its strong 
Markovian assumption it is unsuitable for strongly cou- 
pled superconducting systems. 

In this paper we present a theoretical tool allowing 
to describe dispersive measurements involving nonlinear- 
ities. The tool is developed alongside the example of the 
experimental setup studied in Ref . [2y| . It is based on the 



complex environments approach similar to what is used 
in cavity QED [36j but also in condensed-matter open 
quantum systems |37l [H, [39j] . The idea is to introduce 
the potentially strongly and nonlinearly coupled compo- 
nent of the detector as part of the quantum system and 
only treat the weakly coupled part as an environment. 
In other words, we single out one prominent degree of 
freedom of the detector from the rest and treat it on 
equal footing with the qubit. This "special treatment" 
of one environmental degree of freedom is an essential 
point of this approach because it allows us to describe 
the dynamics of a qubit coupled arbitrarily strong to a 
non-Markovian environment. On the other hand it en- 
hances the dimension of the Hilbert space to be captured. 
This technical complication can be handled using a phase 
space representation of the extra degree of freedom — in 
our example, a harmonic oscillator. 

In Section [IT] we derive the model Hamiltonian moti- 
vated by Ref. [261 ]. For this Hamiltonian we derive in 
Section [Hi] a master equation and present a phase space 
method which enables us to analyze the dynamics of an 
infinite level system. In Section IIVI we demonstrate that 
our method enables to extract informations about the 
measurement process of the qubit, such as dephasing and 
measurement time and also present three different mea- 
surement protocols. 



II. FROM CIRCUIT TO HAMILTONIAN 

We consider a simplified version of the experiment de- 
scribed in Ref. [26| . The circuit consists of a flux qubit 
drawn in the single junction version, the surrounding 
SQUID loop, an ac source, and a shunt resistor, as de- 
picted in Fig.[T] We note here that we later approximate 
the qubit as a two-level system. The qubit used in the ac- 
tual experiment contains three junctions. An analogous 
but less transparent derivation would, after performing 
the two-state approximation, lead to the same model, 
parameterized by the two-state Hamiltonian, the circu- 
lating current, and the mutual inductance, in an iden- 
tical way [l7j . The measurement process is started by 
switching on the ac source and monitoring the amplitude 
and/or phase of the voltage drop across the resistor. 

The SQUID acts as an oscillator whose resonance fre- 
quency depends on the state of the qubit. When the 
measurement is started, the qubit entangles with the res- 
onator an shifts its frequency. The value of the oscillator 
frequency relative to the frequency of the ac driving cur- 
rent determines the amplitude and phase of the voltage 
drop across the resistor. 

The detector (voltmeter) contains an internal resistor. 
This is a dissipative element connecting the quantum me- 
chanical system (qubit + SQUID) to the macroscopic 
observer. The resistor is needed for performing the mea- 
surement, defined as the transfer of quantum information 
encoded in a superposition of states to classical infor- 
mation encoded in the probabilities with which the the 



3 




FIG. 1: Simplified circuit consisting of a qubit with one 
Josephson junction (phase 7, capacitance C q and inductance 
L q ) inductively coupled to a SQUID with two identical junc- 
tions (phases 71,2, capacitance Cs) and inductance Ls- The 
SQUID is driven by an ac bias 7b (t) and the voltage drop is 
measured by a voltmeter with internal resistance R. The to- 
tal flux through the qubit loop is $ 9 and through the SQUID 
is $ s . 



voltmeter shows certain values, e.g. voltage amplitude or 
phase. Note that practically this resistor may be the dis- 
tributed impedance of the coaxial line connected to the 
chip. 

In this section we derive an effective Hamiltonian for 
this system. Our starting point is a set of current con- 
servation equations for the circuit of Fig. [T] 

The total magnetic fluxes through the SQUID and 
the qubit loops can be divided into screening (s) fluxes 
produced through circulating currents and external (x) 



fluxes from outside sources $t 



(x) 



£° and $ 9 



(s) 

$g . Generally, Ampere's law for a system of 



current loops can be represented in matrix form 




which can be inverted 



Ml 



Lq M Sq 

M Sq L s 



Ls -M Sq 
-M Sq L q 




(1) 



(2) 



where Ms q is the mutual inductance and Ms is the deter- 
minant of the inductance matrix M| = L q Ls — M| ? > 0. 
The circulating current through the SQUID loop is given 
by the difference of the currents through the two branches 
Is = [h — h)/2 and the bias current is I = I\ +I2. The 
fluxoid quantization [40| in the two loops reads 



$ 5 

r_ = 71 — 72 = 27T— — mod27r, 

7 = 2n— mod27r, 
$0 



(3) 
(4) 



where $0 — h/2e is the magnetic flux quantum for a 
superconductor. To obtain the equations of motion for 
the phases 7, r± with T± — 71 ± 72 we start from the 



current conservation in each node. 



Ij = I cS sm"fj + VC S 



$0 

= Icssmjj +jj—C s , 
3 G {1,2} 



(5) 



Here we assume that the two junctions have identical crit- 
ical currents. This symmetry, as will be discussed below, 
has the consequence that at zero bias current through the 
SQUID, the qubit will be isolated from its environment. 
In experiment, the two SQUID junction will of course not 
be identical. For an asymmetric SQUID the qubit can be 
protected from environmental noise [4l| by applying an 
appropriate DC bias. 

Using Eqs. and (0 for the circulating current, we 
obtain 



IcS cos — sin — 



if 

Ml \ 



M Sq $ q s) - L q ^ s 



(6) 



Considering the analogy between Josephson junctions 
and inductors we introduce the Josephson inductance 
Ljs = $0/(271" / c s) and rewrite Eq. © using the flux- 
oid quantization ([3]) and ([4]) 



1 r+ . r_ » c s 
_ cos _ sm _ + r _- 



1 

M 



2 (M 8g j-L q T-)+S u (7) 



where the influence of external fields is captured in Si = 
27r/($ M|) (-M Sg $ q x) + L g $P^ . 

For the bias current we have I+V/R = Isit) and from 
Eq. ([5]) we obtain 



Ljs 



■ sin — cos ■ 



f * 



4R 



$0 



I B {t). 



(8) 



For the circulating current in the qubit loop it follows 
from Eq. © 



C q l— + I cq sin 7 



= -^r(Ls<S>i s) -MsA S) 



Ml 



(9) 



Using Lj q = $ /(I cq 2ir), Eqs. ©, and (T5|), this becomes 



Ctfj + - — sin 7 
L>J q 



L s , Ms 



mi" di T ~ + ^ ^ 



—2 



where we defined 

2tt/(M!$ ) (-M Sq <S> ( s x) + L s <Z> q x) ) . 

From Eqs. (J7J), ©, and (fTO]), we observe that T + , the 
phase drop across the SQUID, serves as a pointer: it 
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couples to the qubit degree of freedom 7 and is read out 
by the classical observer, which appears in the classical 
equation of motion ([8|) as a dissipative term. Without 
bias current Is = 0, the classical solution for this degree 
of freedom becomes T + = independent of the internal 
degree of freedom T_ and the qubit. It follows that, 
in the absence of Ib, there is no coupling between the 
quantum mechanical system and the environment, as the 
pointer is decoupled from the observer. 

We start the derivation of the system Hamiltonian sup- 
pressing the dissipative term in Eq. ijHJ). It will be later 
reintroduced in the form of an oscillator bath. Start- 
ing from the equations of motion (JT]), (0), and (fTU)) . for 
r± and 7 we first determine the Lagrangian such that 
dt(djC) = d-yC and dt(df C) — dr ± C. We introduce the 
canonically conjugate momenta p = djC = h^Cqj/e 2 
and P± — <9p ± £ = h 2 C s T±/(2e 2 ) and finally derive the 

Hamiltonian using H = -yp + f + - C. This 

leads to 



H = 



p2 + p2 



C s 



P 

2C a 



e 2 2 r_ r+ 

_ T COS — — COS — — 

h 2 \Ljs 2 2 



M. 



Ml 1 Ml 2 + ^ 



L 



Tq 



■ cos 7 



Lq r 2 

4 

3' 



M% 2 + " 27 " 



-j B (t)r + U.(n) 



Now we proceed to simplify this Hamiltonian using the 
assumptions that Ljs 3> £5, which applies to small 
SQUIDs as the ones used for qubit readout, and that the 
driving strength is small enough to remain in the har- 
monic part of the potential \Ib \ <§C I c s- Using the latter 
assumption, we can expand the potential energies to sec- 
ond order around the minimum and obtain two coupled 
harmonic oscillators (T + and T_) with greatly different 
frequencies. T_ evolves in a much narrower potential 
(oc l/Lg) than that of T + (oc l/Ljg). Therefore we can 
perform an adiabatic approximation and substitute T_ 
through its average position. We obtain the following 
potential for the remaining degree of freedom T + 



U = Un 



1 



cos si 

4Lj S V 2L 



Ml 



Ms. 
' 2L, 



T, - 



I B {t) 



, ^ Mi, M Sq 
C0S ( "121^+7^ 




In the next step, we perform the two-state approxima- 
tion of the qubit along the lines of Ref . [HI, H3| , reducing 
its dynamics to the two lowest energy eigenstates. This 
space is spanned by wave functions centered around two 
values 7 = 7o<7 z (7 either in the left or the right well 
of the potential). While the manipulation of the qubit is 
usually performed at the optimum working point [1 81 ] , the 
readout can and should be performed in quantum nonde- 
molition (QND) measurement i.e. in the pure dephasing 



limit. This reduces the qubit Hamiltonian to eo& z . We 
allow for a significant off-diagonal term oc a x to have 
acted in the past in order to prepare superpositions of 
eigenstates of a z . Physically, this situation is achieved 
by either making one of the qubit junctions tunable, or 
imposing a huge energy bias to the qubit. 

We note here that if, opposed to the case we will dis- 
cuss in the following, the measurement interaction would 
not commute with the qubit Hamiltonian, a full analy- 
sis in terms of quantum measurement theory would be 
required. Similar to Refs. [43|, |44[ the action on the sys- 
tem given by each measurement result would need to be 
determined in order to quantify the information that the 
observer can obtain about the initial state of the qubit, 
as well as the state following the measurement. 

After these approximations the qubit-SQUID Hamil- 
tonian reads 



H S - e(t)a z + -± 
2m 



Pi m(n 2 + A 2 a z , ., 



x 2 - F{t)x , (13) 



where x corresponds to the external degree of freedom 
of the SQUID T + and F(t) — Fosm(ist) originates in 
the ac driving by a classical field. The conversion of the 
parameters to circuit-related quantities can be found in 
Appendix I VII Al Here A is the quadratic frequency shift 
(QFS). 

An important property of this Hamiltonian is the ab- 
sence of the commonly used linear coupling between the 
two-level system and the harmonic oscillator [TtJ ■ In our 
case the qubit couples to the squared coordinate of the 
oscillator, which leads to a qubit dependent change in the 
frequency of the harmonic oscillator instead of the shift 
of the potential minimum. 

Because of the coupling to the driven oscillator the 
qubit energy splitting becomes time-dependent e(t) = 
eo + vI%(t). 

To model the dissipation introduced by the resistor we 
follow the standard Caldeira-Leggett approach [H, |H, 
Efi S3 and include an oscillator bath to our Hamiltonian 



H 



H, 



E 



2m; 



Hsi 



-8{w — u>i) = mhnujQ(uj — w c )/^ 



with J(w) =EiAf- 2milV 
[481 ] where is the Heaviside step function and K = 
the photon loss rate. The cut-off frequency u> c is phys- 
ically motivated by the high frequency filter introduced 
by the capacitors. 



III. METHOD 

Our goal is to analyze the resolution and measurement 
time and investigate the backaction on the qubit. The 
former requires tracing over the qubit and discuss the 
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dynamics of the pointer variable of the detector, the lat- 
ter requires tracing over the detector degrees of freedom. 

It is well established how to do this in principle ex- 
actly [33| when the qubit couples to a Gaussian vari- 
able of the detector (i.e. sum of quadratures of the en- 
vironmental coordinates). A method to map a damped 
harmonic oscillator to bath of uncoupled oscillators with 
a modified spectral density [4!J also exists. In our 
case, due to the quadratic coupling between the qubit 
and the damped oscillator (|13p . any such normal-mode 
transform does not lead to the usual Gaussian model and 
thus cannot use many of the methods developed for the 
spin-boson model. 

There are several approaches to dealing with this chal- 
lenge. As long as the coupling is weak, A <C f2, one can 
still linearize the detector dynamics and make a Gaus- 
sian approximation as it was done in Ref. [28j | . Neverthe- 
less, weak coupling decoherence theory as reviewed e.g. in 
Refs. [3, EiiJ builds on two-point correlators and cannot 
distinguish Gaussian from non-Gaussian environments. 

In this work we describe arbitrarily large couplings be- 
tween qubit and oscillator going beyond the Gaussian ap- 
proximation. The only small parameter we rely upon is 
the decay rate of the oscillator k. This is justified by the 
fact that dispersive measurement only makes sense for 
large oscillator quality factors Q > 1. We treat a com- 
posite quantum system — qubit (g> oscillator — weakly 
coupled to the heat bath represented by the resistor. This 
complex environments approach resembles the methods 
of, e.g. Refs. [H,^. 

We start with the standard master equation for the re- 
duced density operator in Schrodinger picture and Born- 
Markov approximation [32l . l5ll ] , assuming factorized ini- 
tial conditions p(Q) — ps{0) ® Ab(0) 



ips{t) = ^[6s,Ps(t)] 
at m 



1 



(15) 



. dt'Tr B H S B,[Hs B (t,t'),ps(t) <8> /5b (0)] 
(ihY J L 



Ul,H SB Ut and (if 



and T is the time- 



where HsB(t,t' 

Texp^*' dT(H s + H B )/(ih)^ 

ordering operator. In thermal equilibrium there will 
be correlations between the main oscillator and the 
oscillators of the bath, so the initial state is not strictly 
speaking factorized. In the low k limit here, these 
correlations will affect the dynamics only in higher order 
k 2 and can hence be neglected. This is a standard 
assumption in the perturbative treatment of open 
systems where here k is the perturbative parameter, see 
e.g. [13] This approach is valid at finite temperatures 
ksT 3> Hk, for times t 3> l/u> c [SH, [52j, which is the 
limit we will discuss henceforth. We assume unbiased 
noise (Hsb) = 0. The Hamiltonian (flU)) describes 
a driven harmonic oscillator, therefore the Floquct 
modes (see e.g. Refs. [H, [54| for a short review) form 
the appropriate basis in which we express the master 
equation. For a driven harmonic oscillator the Floquet 



modes [54| are given by 

V n (x,t) = tp n (x-£{t)) 



x exp 



mi(t){x-t(t)) 



E n t+ I dt'L(t'] 
'a 



= $ n (x,t)<S 1 



iE n t/h 



(16) 



where E n = Ml(n + 1/2), tp n (x) a number state, £(i) 
is the classical trajectory and L(£, £, t) the classical La- 
grangian of the driven undamped oscillator 



m = 



F sm(vt) 
m(f2 2 — v 2 ) ' 



(17) 



m,it) = ±mt 2 (t)-± m n 2 e(t)+at)F(t)- as) 

We also define the operator A as the annihilation opera- 
tor corresponding to a Floquet mode 



a = A + C(t), 



(19) 



where C(*) = \/mi(i£(t) + so that A$ n {x,t) 

y / rT$„_ 1 (a;, t). After some algebra we obtain 



2mfl 

a(t,t') = <B m ^-^A + c(t') 



(i + it)+^), (20) 
(21) 



Eq. (|2Tj) has been obtained by calculating a(t, i')$ n (x, t), 
where {$„(x, t)} build a complete set of functions at any 
time t. Here one can interpret the sum A + A^ as the 
deviation of x from the classical trajectory £(i). 

Since we are describing a composite quantum mechan- 
ical system, the operators in Eq. (|15[) can be written in 
the qubit a z basis as follows 



Ps 

A = 



_ ( PTT Ptl 

pit pu 
i T o 
o ix 
o e {!,!}, 



a(t,t') 





H sl 



o 



a;(t,t') 



(22) 



-(23) 



where all the matrix elements are operators in the oscil- 
lator Hilbert space 

Asu=±e(()+ (| + ^!±^) i2 - if(() ), (24 ) 



(25) 



and | is the annihilation operator of a Floquet mode 
with frequency S~2-p. j_ = V^ 2 ± A 2 . The functions £(t) 
and £(i) also depend on the frequency of the harmonic 
oscillator, therefore they become 2x2 diagonal matrices. 

As we observed in the previous section, as long as Ib = 
there is no direct coupling between the qubit and the 
oscillator in the second order approximation and thus no 
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coupling to the environment. Therefore, at t = 0, before 
one turns on the ac driving, the harmonic oscillator has 
the frequency f2 independent of the qubit. Therefore the 
initial condition for the density matrix p$ is ps(0) = 

Pqubit <&> PhO- 

We introduce also the annihilation operator Aq for the 
Floquet modes with frequency f2 which relates to A a as 
follows 



A n 



C(t)+ReCo(t)/<r-iImCo(t) 



fa 



(26) 



in in (ihy 



where f a = y/Q a /£l. 



Using the operators introduced above in Eq. ([15]) . we 
obtain 



dt' / dwJ(w)f((B i " (t - t ' ) n(w)+(B" i " (t -*' ) (n(cj) + l)) 
o Jo ^ ' 



[gl (a a + ai) (a a (t,t') + ai(t,t')) p aa , (t) - g a g a , (a a (t, t') + fit(t, t')) Paa ,{t) (a a , + a\, 
[gaga' (a<r + at) pa*>{t) {a a > (t, t 1 ) + a*, (t, t')) - g\, p aa , (t) (a^ (t, t 1 ) + a\, (t, i')) (a<r> +4')) 
(e-^Q-Vniu) + «><*-*'> (nM + !)) }. 



(27) 



where g ff = ^/ h/(2mfl a ) and n(w) is the Bose function. 
We observe that the equations of motion for the four 
components of ps are not coupled to each other. This is 
the consequence of neglecting the tunneling in the qubit 
Hamiltonian Eq. (|13jl . While the two diagonal compo- 
nents fulfill the same equations of motion as in the case 
of the well-known damped harmonic oscillator, each of 
them with a different frequency, the two off-diagonal el- 
ements of the density matrix have a more complicated 
evolution. Specifically, they are not Hermitian and do 
not conserve the norm. This is to be expected, as the 
norm of the off-diagonal elements measures the qubit co- 
herence, which is not conserved during measurement. 

One can handle master equations in an infinite- 
dimensional Hilbert space with the aid of phase-s pac e 
pseudoprobability distribution functions [HI, El, \El 



i.e. 



56, 



which encode any operator with a finite norm 
into a phase-space function. Here, we choose the char- 
acteristic function of the Wigner function x(a, a* ,t) to 
represent the density matrix 



Paa' (*) 



1 



d 2 a Xaa ,(a,a*,t)D(-a), (28) 



where D(—a) — exp(— uA^ + a* Aq) is the displacement 
operator. By replacing this representation of p acr > into 
the master equation (|2T[) we obtain partial differential 
equations for the characteristic functions X ffCT '(a, a*,t). 
Note that here \a) is different from the coherent state as 
it is composed of Floquet modes instead of Fock states, 



|a) = e-H 2 / 2 ]T^=|<I>„(i)). 



(29) 



IV. RESULTS 



Measurement 



In this section we propose three dispersive measure- 
ment protocols, all based on the detection of the oscilla- 
tor momentum, from which the state of the qubit can be 
inferred. 

We start by computing the measured observable of the 
detector, the voltage drop V across the SQUID. In our 
notation, the operator is found as the oscillator momen- 
tum, V — iVbiqubit <£> (a — aJ) , and involves a trace over 
the qubit. Here the momentum p is 2meV/h. Thus, 
we obtain the diagonal characteristic functions Xoa- For 
<7 = a' we obtain from Eqs. (|27|). (|2"gj) . and (|2"oT). a Fokker- 
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Planck equation [5S 



Xaa{a,a* ,t) 



where 



a* (-|-ifi+)+a(-|-ifi-))a a . 
(l + 2n ff )^(a + a*) 2 



Evaluating the integrals in Eq. (|36|) we obtain for the 
probability density of momentum 



(po-C a (t)f 



+ (a + a*)f a (t) 



(38) 

Here, assuming the oscillator initially in a thermal state, 
we have 



(30) C ff (t) 



= iFp (cos(ut) K u + sin(^) (» 2 - » 2 )) 

V2ra~m(n 2 -v 2 ) 



K 2 V 2 + ^2 _ ^2)2 



cos(ut)(fl 2 r — v 2 ) + sin(pt)KP 



+ e-/^ ( n, t )" 3 ^ + tf '-f +"-W+nt (39) 



-e- Kt ^sin(n CT t)- 



2f7 CT (f7 2 - v 2 ) 



and f2± = (±ft 2 + fi 2 )/(2ft). Note that we must ex- 
press the operators in Eq. (fT5|) in terms of Aq corre- 
sponding to frequency Q as the oscillator has initially 
that frequency, in a particular case the thermal state of 
frequency tt. Eq. (|30|) is consistent with the property 
Xaa' (ot) — Xa'a(~ °) originating in the hermiticity of the 
density matrix. We perform the variable transformation 
(a,a*,t) — ► (z,z*,s) defined by means of following dif- 
ferential equations 



d s a = a - i&a+^J + a* 
d s a* = a* (I 



,a = a 
s = t, 



ifl a+ ) + a ( - + ifi ff _ , - 



(32) 

(33) 
(34) 



The solutions of these coupled differential equations will 
depend on some initial conditions i.e. a(s = 0) = z and 
thus we obtain the transformation a — > a(z, z*,s). This 
transformation conveniently removes the partial deriva- 
tives with respect to a and a* in Eq. (|30p and we are left 
with 

dsXaaiz, z*,s) = (a(z,z*,s) +a*(z,z*,s))f (T (t) 

-(1 + 2n ff )^ (a(z, z*,s) + a*(z, z*, s)) 2 , (35) 

which can be directly integrated. After performing the 
transformation back to the initial variables x CTIT (a, a* , t) 
we can calculate the probability density of momentum 
P(Poit) = ^JhmVl/2(8(p — po)) where the qubit initial 
state is q^\ T) + <?ll I) an d 



B a {t) = 



i + 2n(n g ) n a 

2 



(E~ Kt cos(2fi CT t)^2 
4f> 



-hit <y_ 



2fi„ 



l + 2n(0) 



3. 

20? 



1 + — 

O 2 



(40) 



cos(2f7 CT t)- 



(i+S) 



4f> 



4fi„ 



sin(2^)^=-). 



P(p ,t) = nYl 



2n 2 



dk e ifc (S(*) ro -fo) /d 2 « (a, a* , t) 

(36) 
(37) 



°-£{U} 

x [d 2 f3(f3\<e- k ^ A o- A ^D(-a)\(3), 



and f2 ff = y/£l 2 — k 2 /4. One can see that f?(i) evolves 
from B a (0) = l/2+n(0) toB CT (oo) = (l/2+n(f2 (7 ))fi ./ft 
and for f2 = f2 CT (p)cr(t) = C a (t) becomes the momentum 
of the classical damped oscillator with the initial condi- 
tions p(0) = -F Kv/{n 2 -v 2 ) andp(0) = F Q i>/{tt 2 - v 2 ). 
Note that the value B a {po) is independent of the initial 
B a (0) . Therefore the long time value of B is the same 
also for ground and coherent state. 

In the following, when analyzing different types of mea- 
surement protocols, we have to differentiate between dis- 
crimination and measurement time. Measurement time 
is the total time needed to transfer the information from 
the qubit to the observer. In a sample-and-hold protocol, 
one imprints the qubit state into the oscillator, then de- 
couples the two and observes the latter. The time needed 
for the first step is called discrimination time. 



1. Long time, single shot measurement 

In the measurement scheme of Ref. 26:] one needs the 
voltage amplitudes corresponding to the two qubit states. 
For this one must wait until the transients in the momen- 
tum (voltage) oscillations have died out. From the am- 
plitude of momentum one can then determine the state 
of the qubit. 

Following Ref. (26j we define the measurement time as 
the time required to obtain enough information to infer 
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the qubit state 



(n-n) s 



(41) 



where V a is the amplitude of the voltage for the qubit 
in the state \a) and Sy — IksTR is the spectral density 
of the detector output. This is the time needed for dis- 
criminating two-long time amplitudes relative to a noise 
backgroud given by Sy . Therefore, in our notation, 



(42) 



where 



An 



(43) 



and b = ksTC s / 1%. Note that in this type of ampli- 
tude measurement it is advantageous to drive far from 
resonance, since at resonance the amplitudes A a become 
identical for the two qubit states. Off-resonance r m is a 
monotonically falling function of A, i.e. larger coupling 
leads to faster measurement. Close to resonance r m grows 
again for large values of A. 

It is known that, when an harmonic oscillator is driven 
close to resonance, a phase measurement reveals most of 
the information about the oscillator frequency and leads 
to the best resolution and quantum limited measurement. 
Along the lines of Ref. [60] one can suppose that, for 
measurement closest to the quantum limit, the conjugate 
observable to the one being measured should deliver no 
information. In our case, for off-resonant driving and 
amplitude measurement, most of the information about 
the qubit is contained in the amplitude and almost none 
in the phase. 



2. Short time, single shot measurement 

In the measurement protocol of the previous section 
and Ref. [261 ] the desired information is extracted from 
the long-time C a = (p) a (t). The method has the advan- 
tage of being " single shot" , but disadvantages resulting 
from long time coupling to the environment such as de- 
phasing, relaxation and loss of visibility [6l|, [62j . 

In this section we present a different measurement 
protocol. It is based on the short time dynamics il- 
lustrated as follows: for the qubit initially in the state 
1/V2(\ t) + | I)) the probability distribut ion of momen- 
tum is plotted in Fig. [2] (a) and (b). 

In Fig. [2] one can see that the two peaks corresponding 
to the two states of the qubit split already during the 
transient motion of (p)(t), much faster than the tran- 
sient decay time. If the peaks are well enough separated, 
a single measurement of momentum gives the needed in- 
formation about the qubit state, and has the advantage of 
avoiding decoeherence effects resulting from a long time 
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FIG. 2: Probability density of momentum P(po,t) (a), snap- 
shots of it at different times (b) and expectation value of 
momentum for the two different qubit states (c). Here 
hQ,/{k B T) = 2, A/fi = 0.45 , n/Q. = 0.025 and hv/(k B T) = 
1.9 and p is the dimensionless momentum po/vfeT™. 



coupling to the environment. Nevertheless the parame- 
ters we need to reduce the discrimination time also en- 
hance the decoherence rate. 

We define in this case the discrimination time as the 
first time when the two peaks are separated by more than 
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the sum of their widths i.e. 

|Ct( T discr) — (Tdj Bcr ) I 



> 3V2mhn Js T (r disC r) + \ M(r d 



(44) 



A comparison between discrimination and dephasing rate 
will be given in Section HV CI 

Because of the oscillatory nature of C a (t) the problem 
of finding the first root of Eq. (|4~4"|) is not trivial. We solve 
it by semi-quantitatively probing the function | C-p (t) — 
Ci(t)\ - 3V2mHn(y/B^t) + y^JTj) therefore the plot 
ist not very accurate. Nevertheless it gives a good idea 
about the dependence of raiser on A. 
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FIG. 3: Discrimination time as function of the coupling 
strength between qubit and oscillator. Here Ml/(ksT) — 2, 
re/fi = 0.025 and hv/(k B T) = 1.95. 

We observe in Fig.[3]that raiser is a discontinuous func- 
tion of the coupling strength A, such that small adjust- 
ments in the parameters can give important improvement 
of the discrimination time. 

For this type of measurement we are interested in the 
transients of C a {t) and we observe that the difference 
|Cj(t) — Ci(t)\ increases for values of the driving fre- 
quency v close to resonance. For the v far off-resonance 
the splitting of the peaks is increased by stronger driving. 

The discrimination time discussed here is not to be 
confused with the physical measurement time. In par- 
ticular, the discrimination time remains finite even at 
vanishing k and when the off-diagonal elements of the 
full density matrix in qubit space (|22[) still have finite 
norm at this time. The discrimination time is the time 
it takes to imprint the qubit state into the oscillator dy- 
namics. For completing the measurement the oscillator 
itself needs to be observed by the heat bath and, as a 
consequence of that observation, the full density matrix 
will collapse further. 

We note that in this kind of sample-and-hold measure- 
ment, the qubit spends only the discrimination time in 
contact with the environment. Keeping the discrimina- 
tion time short may be of advantage in limiting bit flip 



errors during detection. We do not further describe such 
error processes in this paper. 

As a technical limitation, it should be remembered that 
our theory is based on a Markov approximation for the 
oscillator-bath coupling, hence it is not reliable for dis- 
crimination times lower than the bath correlation time. 



3. Quasi-instantaneous, ensemble measurement 

In this section, we are going to take this idea to the 
next level and analyze a measurement protocol that is 
based on extremely short qubit-detector interaction. In 
Refs. [H, [H, [H[ has been shown that one can measure 
several field observables through infinitesimal-time prob- 
ing of the internal states of the coupled qubit. In this 
section, we apply the same idea to the opposite setting. 
We show that the information about the state of the qubit 
is encoded in the expectation value of the momentum of 
our oscillator at only one point in time, leading to a fast 
weak measurement scheme. 

We rewrite the Hamiltonian (fT4"]) 



H 



m(tfa + - ) + F(t) 



2mQ 



(a + i 



e(t)a z 



HA 2 

' 40 



2 ft (a + at) A; [bi + b 



hioi 



In Schrodinger picture we have 
• 1 



ih 



[H,p], 



which leads, for any observable, to 

1 



d t (o) = (d t o)- 

Setting O — a — a' we obtain: 



([0,H]). 



(45) 



(46) 



(47) 



hA 



nn(a + a^) + —a z (a + a^) (48) 



2mQ 



F(t) 



2mu!i 



-A. 



(bi + bf) 



We assume unbiased noise and the qubit in the pure ini- 
tial state q^\ t) + 9il I) which leads to 



For t = Eq. (J48J) becomes 

d t {a - a*)\ t= Q = ~i{a + a)) t=Q ( 

2\IJ^F{0) 



(49) 



2n 



2mn 

2 



(50) 
(51) 
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If (a + &') t _ 7^ 0, the ensemble measurement of dt(a — 
a^)|t=o is sufficient to determine the state of the qubit. 

In our case the oscillator is initially in a thermal state 
and (a - 



/t=o 



0. Nevertheless, calculating (p) from 
Eq. (|38[) we obtain for the center of the Gaussians corre- 
sponding to the two qubit states C CT of Eq. 



(52) 



(p)(t) = C^ + iC^-C^t))^ 



which is valid for all times. If Cj(i) ^ Cj(t) we have 



( P )(t)-C x ft) 



(53) 



At t = we have like in the exact case Cf(0) = Cj(0) 
independent of system parameters. This is again the 
consequence of the thermal initial state. Therefore one 
cannot infer from (p)(0) the state of the qubit. For in- 
finitesimal r > we have C-f(r) ^ Cj (r). Therefore 
quasi-instantaneous measurement of momentum still de- 
livers the necessary information about the qubit, if the 
measurement is made at a infinitesimally small r > 0. 

At t = also the first derivative of Cf(t) — Ci(t) is 
due to the thermal initial state. A series expansion of 
Eq. (j3"5|) gives the short time result 



<7 T (i)-Ci(t) 



2i/F A 2 
i/ 2 - Q 2 



< 2 +C(i 3 ). 



(54) 



This gives a criterion for T,n scr , independent of n, similar 
to observations of previous section, i.e. for short dis- 
crimination times we need large A and strong, close to 
resonance driving. 

Moreover, it is sufficient to measure the expectation 
value of momentum, and not the first time derivative. 
The reason for this is the oscillator evolution, mediated 
by the interaction with the qubit, into a state with fi- 
nite expectation value (a + af), in other words the sys- 
tem is automatically creating its own measurement favor- 
able "initial" condition. This is visible in Eq. where 
the part of the signal proportional \q^\ 2 increases like t 2 
which, for short times is slower than t, as it would be 
in the case where the favorable initial condition already 
exists. 

This method leads to shorter discrimination times than 
the protocols presented in section IIV A II and IIV A 21 
which arc independent of A and v. Again, the read out 
of the oscillator in the end will be a separate issue and 
ultimately take a time oc k _1 . 

On the other hand, in order to establish the expecta- 
tion value with sufficient accuracy, this scheme requires 
a large ensemble average. According to the central limit 
theorem the uncertainty of the ensemble measurement is 



Ay 

(y) 



1 Ap 

n(p) 



(55) 



where y is the ensemble averaged value of the measured 
momentum and N the number of measurements. For a 



given precision we have N oc (Ap/(p)) , therefore the 
number of measurements necessary to reach a given pre- 
cision depends on A and time t. We have 



Ay 2 



(E. 6{U } l««r| 3 C«r(t) 



Eq. (|55|) shows that we need Cj ^ Cj. in order to de- 
termine the state of the qubit. At the same time, the 
number of measurements N necessary for high precision 
measurement of momentum is significantly reduced when 
Cf = Cj (the two Gaussian distributions overlap com- 
pletely). This reflects the tradeoff between the number 
of measurements and the signal strength Cf — Cj which 
provides the information about the qubit. 



B. Back-action on the qubit 

In order to complete the study of the measurement 
protocols presented in the previous section, we need in- 
sight into the measurement bakaction on the qubit. Since 
we are studying the QND Hamiltonian (fl~3| . the qubit 
decoherence consists only of dephasing. We start with 
the qubit in the initial pure state (| t) + I J.)) A/2 and 
study the decay of the off-diagonal elements of the qubit 
density matrix. Such a superposition can be created by 
e.g. rapidly switching the tunnel matrix element from a 
large value to zero [661 ] . or by ramping up the energy 
bias from zero to a large value. We compute the qubit 
coherence 



C{t) = Tr(a x ® lps(t)) =2ReTrp n (t) 

= 2Re J dxj dpW n (x,p,t) 

= 2Re JdxJdpe ix0 e ip0 W n (x,p,t) 
= 8^Re XT i(0,0,t), 



(57) 



where W aa i is the Wigner function 



W atT '{x ,Po) = / dy(x +y\paa>e 2iypo \xo - y)- 



(58) 



We extract the dephasing time from the strictly expo- 
nential long-time tail of C(t). 

We rewrite the master equation (f2"T|) for a ^ a' us- 
ing Eqs. (|2"6")) and l|2"8j) and obtain a partial differential 
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equation for the characteristic function xn 



equations of motion 



A E {t) = Bs\n{i>t){R{t) - S{t))+T{t) 



(a(ki + iQ) + a*ki+ B sin(^i)) d a 
+ (a* (k 2 - ifl) + ak 2 — B sin(t4)) d a * 

~ i §j( 5 a- a a*) a + (« + «*)/u(*)+^(*) 



+ p{a + a*) 2 ) xt J. (a, a*,*) 



(59) 



The coefficients can be found in Appendix IVII Bl 
Eq. is a generalized Fokker-Planck equation where 
the total norm is not conserved, i.e. Jd 2 ax^i{a,a* ,t) is 
not a constant of motion. 



1. Solution of the generalized Fokker-Planck equation 



Generalized Fokker-Planck equations (f59"|) cannot in 
general be solved analytically with the established tools 
[53 |. In our case, we are not interested in a fully gen- 
eral solution of the differential equation, but in the ini- 
tial value problem where the xn(a, a*,0) is a Gaussian 
function, which covers thermal and coherent states of the 
oscillator. In this case one can show that Xn( a i °* ' i ^) re ~ 
mains a Gaussian at all time. This is the consequence of 
the QND Hamiltonian (flU|) . We make the ansatz 



R(t) 
S(t) 
M(t) 
N(t) 

Q(t) 



i§!(M( i) + WW - W -ffi«)(6 1 , 

(fci + ifl)R(t) + k 2 S(t) - Bsm{vt)(2M{t) - Q{t)) 
iA 2 



n 



■(R(t) - S(t))(Q(t) - 2M(i)) + f n (t) (62) 



= (k 2 - iQ)S(t) + kiR(t) + Bsm(vt)(2N(t) - Q(t)) 

- I A-( S (t)-R(t))(Q(t)-2N(t)) + f n (t) (63) 

= 2(ki + iO)M (t) + k 2 Q(t) 

+ 1 -^(Q(t)-2M(t)f-p, 

= 2(k 2 -ifl)N(t) + k 1 Q(t) 

+ 1 -^(Q(t)-2N(t)f-p, 

= (fci + k 2 )Q(t) + 2k x M{t) + 2k 2 N{t) 
iA 2 



(64) 



(65) 







-(Q(t)-2M(t))(Q(t)-2N(t))-2p, 



(66) 



where A(t) = <£ AEi ' t \ This system of equations can be 
solved numerically, for example using a Runge-Kutta al- 
gorithm. 

Ref. ijTij gives a elaborate analysis of the various de- 
phasing mechanisms in the case without driving and the 
parameter regimes where they come to play. There the 
weak qubit-oscillator coupling (WQOC) regime is asso- 
ciated to a phase Purcell effect [68| where the dephasing 
rate 1/r^ oc 1/k. Beyond the weak coupling, Ref. [67| 
explores a strong dispersive coupling regime with fun- 
damentally different origin where the dephasing rate is 
proportional to k. 

In the following we want to apply and extend this re- 
sults to the case of actual measurement, i.e. when the 
oscillator is driven in order to measure its frequency and 
from this information, to infere the state of the qubit. 



Xu (a, <**,*) = A(t) exp ( - M (t)a 2 - N(t)a* 2 
- Q{t)aa* + R(t)a + S(t)a*), 



(60) 



and obtain for the time-dependent parameters of the 
Gaussian a closed system of nonlinear differential equa- 
tions of the first order, thus proving that our ansatz is 
correct and complete if the initial characteristic function 
is a Gaussian. 



2. Qubit dephasing 

We start by studying the dependence of the qubit de- 
phasing on the parameters of the oscillator driving field. 

In Fig. [4] we observe that the dependence of the de- 
phasing rate l/r^ on Fq is quadratic. For values of k 
belonging to strong and weak coupling regime at Fq = 
we obtain the same driving contribution to the dephasing 
rate, proportional to kFq, see the inset of Fig. |4j Here 
only the contribution of driving is shown. We have sub- 
stracted from each curve the initial value of at Fq = 0. 
We observe that the decoherence rate must be of the form 



1 



1 



ct. • F q 2 k, 



F o =0 



We assume the oscillator initialy in a thermal state 
Xu(a,a*,0) = (l/47r)exp(-(l/2 + n(tt))H 2 ) and ob- 
tain for the parameters of the Gaussian ansatz following for both the weak and strong coupling regime. 



(67) 
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FIG. 4: Dephasing rate dependence on driving: dependence 
on A for different driving strengths Fo (k/Q = 1CF 4 and v — 
2fl). Top inset: dependence of the decoherence rate on Fo for 
different values of re (A/fi = 5- 10~ 2 and v = 2f2). Bottom in- 
set: dependence of the decoherence rate on driving frequency 
v for different vales of re (A/fi = 0.5). Here ftfi/fcsT = 2 and 
Fo is the dimensionless force FqH/ (fcsTvmfcsT 5 )- 



This was to be expected since the qubit couples to the 
squared coordinate which (at least in the classical case) is 
proportional to Fq . In both regimes, the driving leads to 
a contribution to the dephasing rate that is proportional 
to k because the driving leads to classical motion relative 
to the heat bath, which is fixed in the cc-coordinate space. 
This motion enhances the effect of the bath the stronger 
the friction coefficient K is. Consequently, even if in the 
undriven case the dephasing rate scales as 1/k, strong 
driving can in principle cross it over to a decay rate cx n. 
This cross-over from 1/k to k inside the WQOC regime 
happens at either very strong driving or when the driving 
v frequency approaches one of the system resonances fi CT . 

The dependence on the driving frequency has also been 
analyzed in Fig. |4j Here we observe two peaks at f2j 
and f2j. At v — f2 the classical driven and undamped 
trajectory £(t) diverges. In terms of the calculation this 
means that the Floquet modes are not well-defined when 
the driving frequency is at resonance with the harmonic 
oscillator — we have a continuum instead. Physically 
this means that at t = our oscillator has the frequency 
f2 because it has not yet "seen" the qubit, and we are 
driving it at resonance, and by amplifying the oscillations 
of (x) which is subject to noise we amplify the noise seen 
by the qubit. The dephasing rate is also expected to 
diverge. The peaks at J7|and fij, show the same effect 
after the qubit and the oscillator become entangled. The 
dephasing rate drops again for large driving frequencies 
to the value obtained in the case without driving. 



C. Comparison of dephasing and measurement 
times 

In this section we analyze the measurement times nec- 
essary for the measurement protocols described in sec- 
tions IIV A II and IIV A 21 and compare them with the de- 
phasing times of the qubit obtained for the same param- 
eters. 



1. Long time, single shot measurement 

For the long time measurement protocol (section 
IIV A II) we observe that l/r m cx F 2 A 4 + 0(A 8 ). 




io" 3 io" 2 A/a io" 1 



FIG. 5: Comparison of dephasing / measurement rate. Open 
symbols: l/(7"^,fi) and filled symbols: l/(r m fi). Here v = 2fi, 
re = 10 _4 fi and HQ/ '(ksT) — 2. Inset: dependence on re, 
Fo = 200. The crossings, conflicting with the quantum limit 
[l3 |. are signaling the limits of the Born approximation, as 
described in text. 

Comparing l/r^ and l/r m we find that the measure- 
ment time depends more stronlgy on the driving strength 
Fo than the dephasing time. 

As one can see for the parameters of Fig. [5l in the 
WQOC regime the measurement time is longer than the 
dephasig time. Their difference decreases as we increase 
A due to the onset of the strong coupling plateau in the 
dephasing rate, approaching the quantum limit where 
the measurement time becomes comparable to the de- 
phasing time. Note that, for superstrong coupling cither 
between qubit and oscillator or between oscillator and 
bath, corrections of the order (re/fi^) 2 of the dephasing 
rate gain importance. These corrections are not treated 
in our Born approximation. Therefore the regions where 
the dephasing rate becomes lower than the measurement 
rate, in violation with the quantum limitation of Ref . [l4j , 
should be regarded as a limitation of our approximation. 

The inset in Fig. [5] shows the dephasing and measure- 
ment times as function of re. Again we observe improve- 
ment of the ratio of measurement and dephasing time as 
we increase A. On the other hand, if the tunning of re 
should be easier to achieve experimentally, we also see 
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that, at given A one can make use of the phase Purcell 
effect, which reduces the dephasing rate as 1/k while the 
measurement rate increases like k. This goes along the 
lines of Ref. [67| where it has been shown that strong k 
implies WQOC, i.e. phase Purcell effect. 



2. Short time, single shot measurement 

As already mentioned, for the short time, single shot 
measurement strong, close to resonance driving is needed 
for the rapid separation of the peaks. While the discrim- 
ination time is not very sensitive to the change of k, we 
observe in Fig. [5] that one needs relatively strong cou- 
pling (A /CI E (0.03,0.1)) for the discrimination time to 
become shorter than the decoherence time. The picture 
of the dephasing rate is also qualitatively different from 
the case without driving or with far off-resonant driving, 
since for A/f2 ss 0.4, fif becomes resonant with the driv- 
ing frequency. In this region our numerical calculation 
also becomes unstable. Nevertheless, as one can see in 
Fig. the dephasing rate is proportional to k. Thus, by 
reducing the damping of the oscillator, one can extend 
the domain of values of A where the measurement can 
be performed. We also observed that by further reducing 
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FIG. 6: Comparison of dephasing / discrimination rates. 
Open symbols: 1/(t$Q) and filled symbols: l/(r<iiscrf2)- Here 
v = l.lfi ,F = 200 and h{l/(k B T) = 2. Inset: dependence 
on k for A/n = 0.1. 



A the discrimination rate suddenly drops to zero, i.e. for 
too small A the two peaks in Fig. [2] will never be well 
enough separated to allow a single shot measurement. 

In this protocol what we call "discrimination time" is 
actually the time when the system becomes measurable, 
i.e. one can in principle extract from a single measure- 
ment the needed information about the qubit (and con- 
sequently collapse the wave function) . We do not further 
describe this collapse here. 



V. CONCLUSION 

We have presented a phase space theory of the mea- 
surement and measurement backaction on a qubit cou- 
pled to a dispersive detector. 

We have studied the qubit coupled to an complex en- 
vironment (weakly damped harmonic oscillator) with a 
quadratic coupling, which does not have to be weak. We 
solved the problem by considering the prominent degree 
of freedom of the environment, i.e. the main oscillator 
as part of the quantum mechanical system and explicitly 
solving its dynamics. Only at the end of the calculation 
we traced over this last degree of freedom of the environ- 
ment in order to obtain the qubit dynamics. 

We presented three measurement protocols and com- 
pared the measurement and decoherence times. The pro- 
tocol of section IIV A II requires long measurement time, 
such that the measurement can only be preformed in 
strong coupling regime, with far off-resonant driving. 
The protocol of section HV A 21 has the advantage of short 
discrimination times compared with the dephasing time, 
requires strong qubit-oscillator coupling and also close to 
resonance driving. Both this protocols can be performed 
as single shot measurement, and thus may be useful as 
a readout method for the scalable architecture with long 
range coupling using superconducting flux qubits [69| . 
The quasi-instantaneous measurement protocol of section 
IIV A 31 has the advantage of the shortest possible discrim- 
ination time and no restriction for the qubit-oscillator 
coupling, with the drawback that one needs to repeat 
the measurement a large number of times to obtain the 
momentum expectation value. 

We expect our results, with minor adaptations, to be 
applicable to various cavity systems, e.g. quantum dot 
or atom-based quantum optical schemes (70, [zH- The 
dispersive coupling of Hamiltonian (| 13[) could have im- 
plications for the generation of squeezed states, quantum 
memory in the frame of quantum information processing, 
measurement and post-selection of the number states of 
the cavity. 
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VII. APPENDIX 

A. Parameter conversion 

Here we present the parameter conversion from the 
actual circuit to our model Hamiltonian. From the 
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Description 


Symbol 


Circuit 


HVom i on r*\r 


Q 


J * 

V LjsCs 


Mass 


m 


h' z Cs 

1 2 * 1 


Qubit coupling 


A 


/ ^f, 
/ "Jl 

V L JSC S 


Driving field 


F(t) 




Qubit driving 


<t) 


e + vl%(t) 


Momentum 


P+ 


ti"C s -h C S HV 
2,2 1 + - *e 


Position 


X 


r+ 


Damping constant 




1/(2RC S ) 



TABLE I: Parameter conversion 



relation /( 7 ) = /( 7o <t 2 ) = (/( 7o ) + /(- 7 o)) /2 
(/(7o) - /(-7o)) /2 we have: 



B. Parameters for the generalized Fokker-Planck 
equation 



h,2 
P 
B 



l + 2nj)n ; 



K (l + 2n T )fi T 
"2 4f7 

(fl T (l + 2n T ) + + 2nJ) - 

2LF A 2 



— v*- 



1 / / a + 6 7o 
/i = t: cos 

2 V V 2 

„ 1 / ( a + 6 7 o 



cos 



COS 



a - fryo 
2 

a - 6 7o 



2L n 



17 — r^ 2 



Sf 2 



a + 6 7o 



a - 6 7o 



/u(*) 



2^(J-^) 2 (cos(2^)-l)-2-eft) 



sm(vt) | 



A 2 



E 



CT (l+2n CT ) 



^A 2 (l + 2n„) 
(^ 2 -f7 2 )(i/ 2 -ft 2 ) 



ft 2 - z/ 



■ cos(i4) 
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